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高 维 抛物 型 方程 的 一 个 高 精度 恒 稳 定 的 交替 方向 格式 * 


Jk Æl, ERR? 
(1- 福建 生物 工程 职业 技术 学 院 基 础 部 ， 福 州 350002; 2- 华侨 大 学 数学 科学 学 院 ， 福 建 泉州 362021) 
摘 要 : 本 文 对 高 维 抛物 型 方程 构造 了 一 个 高 精度 恒 稳定 的 交替 方向 格式 ， 当 2 < p < 4 时 ， 格 式 的 局 
部 截断 误差 阶 可 达到 O (T? 十 ht)， 当 p = 3 时 ， 将 Richardson 外 推 法 应 用 于 本 文 格式 ， 得 到 
T O(r? + hE) 阶 精度 的 近似 解 . 最 后 通过 数值 实例 验证 了 我 们 对 所 得 格式 所 作 的 理论 分 析 是 正 


确 的 . 
关键 词 : 高 维 抛物 型 方程 ， 交 替 方 向 格式 ， 高 精度 ; 恒 稳定 
分 类 号 : AMS(2000) 35A05; 65M12; 92C37 中 图 分 类 号 : O241.8 文献 标识 码 : A 
1 引言 


在 研究 热传导 过 程 、 气 体 扩 散 现 象 和 电磁 场 的 传播 等 问题 时 ， 常 遇 到 已 维 抛物 型 方程 . 本 
文 研究 如 下 抛物 型 方程 
ðu P. 8u 


a EN O«czi«L, l= 1 2 ,p, t0 (1) 
l=1 "i 


的 第 一 边 值 问题 . 

用 差分 方法 求解 上 述 方程 ， 目 前 已 经 有 了 一 些 较 好 的 格式 [1， 但 有 的 是 条 件 稳定 ， 有 的 是 
精度 不 够 高 . 求解 高 维 抛物 型 问题 比较 理想 的 是 用 交替 方向 法 . 它 是 将 高 维 问题 转化 为 一 系 
列 的 一 维 问题 ， 通 过 一 系列 的 三 对 角 方 程 组 的 求解 来 得 到 高 维 问题 的 数值 解 ， 此 方法 最 早 是 
由 Peaceman 和 Rachford 提出 的 外. 它 每 一 步 仅 解 一 个 三 对 角 线 方程 组 ， 且 保证 格式 的 绝对 稳 
定性 ， 将 显 格式 的 计算 简便 性 与 隐 格 式 的 绝对 稳定 性 结合 起 来 . 但 P-R 格 式 仅 适用 于 求解 二 
维 抛物 型 方程 ， 不 能 将 它 推 广 到 三 维 以 上 的 抛物 型 方程 . 现 有 的 三 维 ADI 格 式 大 多 数 精度 都 
较 低 ， 例 如 Douglas-Gunn 格式 精度 为 O(r2 + h?2). 为 了 提高 精度 ， 本 文 将 Douglas-Gunn 格式 
加 以 改造 ， 使 之 保持 了 原 有 ADI 格 式 的 计算 简便 性 与 绝对 稳定 性 ， 同 时 将 局 部 截断 误差 阶 提 
高 到 O(r2 + ht). M42 «p < 4 时 ， 本 文 给 出 的 格式 是 绝对 稳定 的 . 另外 ， 为 了 进一步 提高 精 
度 ， 我 们 将 Richardson 外 推 法 应 用 于 所 得 的 格式 ， 得 到 了 有 具有 O? + h) 阶 精 度 的 近似 解 . 文 
末 的 数值 例子 表明 我 们 所 作 的 理论 分 析 是 正确 的 . 


2 差分 格式 的 构造 


设 7 为 时 间 步 长 ，Azi (1 = 1,2,… ,p) 分 别 为 z1(1 = 1,2,… ,p) 方 向 的 空间 步 长 .为 简 
便 计 ， 本 文 取 等 步 长 Ax = E = h(M AE), 方程 (1) 的 解 函数 为 u(z1, zz, … Ept) 
id u(jAh, j2h,: 23  jph,nT) = u(ji,j2,: Eo , jp, n)- 
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M E T 二 oo ëa 


用 如 下 的 ADI 格 式 逼 近 方 程 (1)， 有 


p pP 
E -7G - dnt ed a Y 8h8 we, 
e E (2) 
1 


[I — (5 — 33-)82.](w nt Lam) att; oam, s = 2,3,... ,p, 


其 中 j. j, 表示 在 节点 (jih, jah,- ,jph, i Crus 为 简化 记号 ， 在 以 下 的 
os RSS qn net dn, MU ERU in -a 2,( = 1,2,- …,p) 分 别 为 关于 zi1(l = 
… p) 的 二 阶 中心 差 商 ,，r = z 为 网 格 比 . A (2) 的 截断 误差 阶 ， 先 消去 (2) 式 中 
ee ig 得 到 如 下 格式 
1 1 z 
IIt-- gos i)4 ,| tl j=e( 5 + UE. 2- 6282, )u ; (3) 
格式 (2) 和 格式 (3) 是 等 价 的 ， 二 者 具有 相同 的 截断 误差 阶 与 稳定 性 . 
将 (3) 式 中 各 个 节点 上 的 飞 以 其 在 节点 (jih, jah,- ,Jpj nT) 处 展开 的 Taylor 级 数 代 入 ， 并 
注意 到 方程 (1)， 可 得 


Ou | TrOwu 1 1,2100u rB u u 
atza -3(" 5) * 4 (22 agog - o) 
ick 
alo a IPE ALEA 
Otüx? —, OtOr? Ox? 
<k 


_ Ou hu Ou 


= ót 十 1288 08 m 2». 6zg ar DE BG =, Oko 3:153) 十 O(A$). (4) 


于 是 由 上 式 可 知 格式 (3) 或 与 之 等 价 的 (2) RE Ras O(r? 十 2 


3 稳定 性 分 析 
根据 差分 格式 稳定 性 分 析 的 Fourier 方法 ， 令 
u? vus =p Peili +japat+ipbp)h j= VTI, (5) 
且 记 
6 一 —4s2 uj, jay jp? l=1,2,--- ,p, (6) 
其 中 


qur 
m EC 1 = sx = P = 
sin 2M € [0,1], 121,2, ,p, qu=0,1, ,M-—1 


将 (5) ARA 全 式 ， 经 计算 整理 ， 并 利用 (6) 式 的 关系 可 得 格式 (2) 传播 因子 为 


p 
II [14r $-a3)4&]-— 人 37 E sas A 


psi 


p= 7 , (7) 
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其 中 
QUU . Qk7 
Sa = Sino t, So = Sin F717 € [0,1], Lk-—L2,-,p qu5gdk—0,1,---,M — 1. 
要 使 两 层 差 分 格式 (2) 稳定 ， 须 满足 Von Neumann 条 件 |p| < 1， 即 只 须 证 明 
p p 
2 4 2,2 
2r 22 Sa 3T 2 Sq, Sqr 
lk « 1, (8) 


Oi 


p 
Ties - 324] 


成 立 就 行 。 
不 等 式 (8) 的 左 端 只 当 p < 4 时 成 立 ( 左 端 不 等 式 与 7 无 关 )， 实 际 上 ， 因 (8) 式 左 端 中 分 母 恒 


为 正 ， 故 只 要 分 子 为 正 即 成 立 . 


当 p = 2 时 ， 有 
9 2 2 4 2 2 >2 2 2 4 2 >0 
(sa 十 Saz) F 3?« $9, 二 (sq + 52.) p. 3?n 
当 p 二 3 时 ， 有 
4 
2 2 2 2 aa 2 .2 2.5.2 
2(sz, F 302 二 Sg) dii 3 (Sa Sg T Sq2Sq3 + 82:50.) 
4 
2 2 2 2 2 2 
2 2(s?， + Sga 十 Sg) 一 AGA + $$, + Sg) z 0. 
当 p = 二 4 时 ， 有 


4 
2 2 2 2 522 2 2 2 .2 2 .2 2 .2 2 .2 
2(s2, + Sga T Sga + Sg 3 (sq $0, T Sq Saa 十 S4, 5g4 + 99,994 十 95,99, + Sas 52.) 


2 2 
将 3 5 (sz — Sa) 20. 


pr 
不 等 式 (8) 的 右 端 ， 右 端 不 等 式 与 rx 有 关 . 
1 2 2 2 1 1 2 2 .2 r 
t= 了 (sa ME Sga) + 16r (5 一 =) $9189; 十 g Sq Sa 
= 1 1 2 1 1 2 4 2,2 2 -0 
= uu vU ERIT Z O 
Mp 一 3 时 ， 有 
1 2 2 2 211 IX 6n 2 .2 2 .2 
1— 了 (su 十 $$ S2.) t 16r (5 = z) (52,52, E Sgi Sga 十 52,54.) 
1 1 3 4 2 2 ,2 2 2 
+64r? (7 d zx) 52, 82,52; * 37 (sa $0 十 Sg Sgs F $9; 8g; 
1 2 1 2 1 2 2/(,2 .2 2.2 2 .2 2,2 ,2 
z (1 ili 354) (1 7 355) (1 Is 352) +4r (sa, Sa + Sq Sq T 804894 一 51.542521) 


2 
+ (se? 4 37) 82, 52,52, 2 0, 
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1 14? 
2 2 22 
1-335 (ie (5- m) 25 6 
2 .2 .2 1 My 2,4 : 2 .2 
464r? (3 -u) > Sq: Sa; Sqn + 296r ‘(3 一 =z) Is. T 3" 2. 50,53; 
4 4 8 4 4 4 
c 2 22,2 2 2 22 2,22 
a II (1 "3 7Sa ) +r(3 » di dj dk z 115) 十 47 ( 5 Sqi Sa; ET 2: $9,925 Sar 
Minis i272 vu T: viie 


4 
+5 JJs2)+r(8 5 82,52 am - STIS) es TE > 0. 


由 于 格式 (2) 和 格式 (3) 是 等 价 的 ， 因 此 我 们 得 到 如 下 定理 . 
定理 1 当 2<p< 4 时，ADI 格 式 (2) 是 绝对 稳定 的 . 


4 ”外 推算 法 


为 了 进一步 提高 格式 (2) 的 数值 解 的 精度 ， 可 以 应 用 Richardson 外 推 法 于 ADI 格 式 (2). 例 
如 ， 根 据 第 2 节 的 论述 ， 当 p = 3 时 ， 我 们 用 w? pg 表示 格式 (2) 按 空间 步 长 h， 时 间 步 
长 7， 计 算 所 得 的 在 第 nv 层 的 (j1h, jah, jah) 处 的 数值 解 作为 微分 方程 (1) 的 精确 解 ula, jo, a, n) 
的 近似 .由 (4) 式 可 得 


u(j1, j2, ja; n) TE V5, jaja (h, T) 


= T?P(zi, £2, £3, t) + hQ(zi,22,23,t) + O(T? + 7?h? + rh^ + hê), (9) 
其 中 
lrO, ôu Ow 0*u Ou 
P(z1, 22, 23,1) = 4 don i s ü 8353) B o8] 
tu tu tu tu 
Qs, 22, 23,1) = za x E "nd Ox i az?6z3 0z20z2 5:522) 
-x23 (Du $3) + 人 06u $ Su 06u 
Ozi0x2 rðr? + E 
dfu " Ou à Ou 
Oz20z$ Ox20z4 pd 
因此 有 


NP hr 
u(j1. j2, ja; n) z UIS, 2ja,2ja (5; 7) 


- i ? P(zi, 22, 23, t) + ; Q(z1, £2, 23,t) + O(r? + 7?h? + Tht + hô). (10) 
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ey 


EE: hc 
u(j1, j2, ja, n) 一 15 E [16ufs. 22,2js (5; 7) — U}, ja is Ub 7) = O(r? + T? h? + Tht + hê). 


因此 ， 若 用 


zl 6u$ ,272， as (7 T) - V5, ja ja Us 7) 


EX ul jz ja, n) 的 近似 ， 同 时 注意 到 {O(r2j2),O(rj2)} € max{O(7”),O(hs)}， 其 误差 阶 
可 提高 至 O(r3 + hô). 


5 ”数值 例子 
在 区 域 DD : {0 < zi,za,za < T,t > 0} 上 考虑 如 下 初 边 值 问题 
u(z1, o T3, ð = singısinzsinz3, 
u(0, £2, £3, t) = U(T, 22, £3, t) = 0, 
u(x1,0, 23, t) = u(z1, T, 23,0) = 0, 
u(21,22,0,£) = u(z1, £2, T, t) = 0. 
下 面 我 们 给 出 本 文 格式 (2) 24 p = 3 时 的 数值 解 与 精确 解 


—3t 


u(271,22,23,0) = e " sin z1 sin x2 sin za 


和 文献 4 "P Douglas-Gunn 格式 的 数值 解 的 比较 ， 同 时 给 出 了 当 p = 3 时 格式 (2) W9 Richardson 
外 推算 法 的 数值 解 . 这 里 取 Azl = Ax; = Art = h = Eo T= rh, 7 = 者 和 1. 计算 
$jn = 200， 其 数值 结果 如 表 1. 


表 1: 各 种 格式 计算 结果 比较 





























(ji, j2, ja) 格式 (2) 文献 四 格式 外 推算 法 
7.041774e-008 7.043219e-008 7.308086e-008 7.041769e-008 








3.352972e-006 3.353660e-006 3.479778e-006 3.352969e-006 


3.181644e-011 | 3.425410e-011 | 3.179807e-011 


5.169982e-011 5.172921e-011 5.569250e-011 5.169934e-011 
8.485387e-011 8.490212e-011 9.140699e-011 8.485309e-011 







































从 以 上 结果 可 以 看 出 ， 本 文 格式 的 数值 解 与 精确 解 有 较 好 的 吻合 ， 它 较 文献 [和 的 Douglas- 
Gunn 格 式 的 解 至 少 精确 两 位 有 效 数字 . 本 文 格式 当 p = 3 时 ， 外 推 一 次 后 所 得 的 解 的 精度 又 
有 明显 的 提高 ， 这 表明 我 们 所 作 的 稳定 性 分 析 和 精度 分 析 是 正确 的 . 
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An ADI Scheme with High Accuracy and Absolute Stability for 
Solving High Dimensional Parabolic Equations 


ZHANG Xing!, SHAN Shuang-rong? 


(1- Department of Fundamental Courses, Fujian Vocational College of Bioengineering, Fuzhou 350002; 


2- School of Mathematical Sciences, Huaqiao University, Quanzhou, Fujian 362021) 


Abstract: In this paper, an ADI scheme with high accuracy and absolute stability is constructed 
for solving high dimensional parabolic equations. As p varies from 2 to 4, the scheme is absolutely 
stable and the truncation error can reach the order of O(7? + h^) at the maximum. When p is 3, the 
Richardson's extrapolation method is successfully applied to the scheme and the approximate solution 
with accuracy O(r? -- h$) is attained through an extrapolation. The correctness of the proposed scheme 


is verified by numerical examples. 
Keywords: high dimensional parabolic equation; alternating direction implicit scheme; high accuracy; 


absolutely stable 
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